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Abstract 

Critical scaling and universality in short-time dynamics for spin models on a two- 
dimensional triangular lattice are investigated by using Monte Carlo simulation. Em- 
phasis is placed on the dynamic evolution from fully ordered initial states to show that 
universal scaling exists already in the short-time regime in form of power-law behavior 
of the magnetization and Binder cumulant. The results measured for the dynamic and 
static critical exponents, 0, z, (3 and confirm explicitly that the Potts models on the 
triangular lattice and square lattice belong to the same universality class. Our critical 
scaling analysis strongly suggests that the simulation for the dynamic relaxation can 
be used to determine numerically the universality. 
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1 Introduction 



The study of phase transitions and critical phenomena has been an attractive and impor- 
tant topic in statistical physics for a long time [0]. As is well known there appears critical 
universality due to the infinite spatial and time correlation lengths at the continuous phase 
transition points, and the universal behavior of a critical system is characterized by a number 
of critical exponents such that the models in the same universality class have the same values 
of the critical exponents. Therefore, an estimation of critical exponents to determine the 
universality class for different statistical models is an interesting challenge. However such 
research is difficult to be carried out analytically since only few statistical models can be 
exactly solved [0]. Therefore, numerical simulations supply a powerful method to estimate 
critical exponents which are measured by generating equilibrium states and averaging over 
independent configurations in Monte Carlo (MC) studies. For example, the dynamical ex- 
ponent z can be measured from the exponential decay of time correlation for finite systems 
in the long-time regime. Unfortunately, MC simulations near the critical point in the equi- 
librium suffer from critical slowing down (CSD) because the auto-correlation length r 
divergs as r ~ L^, where L denotes the linear size of a lattice and z > 2. 

In the last decade, the exploration of critical phenomena has been greatly broadened and 
much progress has been made in critical dynamics. An important discovery is that there 
exist critical scaling and universality for systems far from equilibrium Traditionally, it 
was believed that universal behavior exists only in the equilibrium or in the long-time regime 
of dynamical evolutions. However, recent researches in critical dynamics for many statistical 
models have shown that universal scaling behavior also emerges within the macroscopic short- 
time regime of dynamic processes after a microscopic time scale tmic 0) |§] • The investigation 
of the short-time critical dynamics (STCD) not only exhibits the existence of universal 
dynamic scaling behavior within the short-time regime, but also supplies a very efficient 
method to determine the critical exponents. 

In Ref.|^, the recent progress in short-time dynamical studies has been reviewed, where 
several simulation methods and the corresponding numerical results have been presented. 
Up to now, the values of critical exponents for a variety of statistical models, such as the 
two-dimensional (2d) Potts and 0^ models on a square lattice 0; |§) 0]; the 2d quantum XY 



and the random-bond Potts models on the square lattices ||10| , |TI|, the 2d FFXY model and 



2+l)d SU(2) lattice-gauge models 0, |T3l, have been calculated precisely. But only very 
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few attention has been concentrated on the 2d triangular lattices |]T^, [15 . 

By the STCD method we can estimate not only the dynamic critical exponent 2;, but also 
the static critical exponents 6', (3 and u. More important, the results for both the dynamic 
critical exponent z and static critical exponents, (3 and u, are consistent with those obtained 
by the traditional MC simulations performed in equilibrium. Furthermore, similar to the 
measurements of critical exponents, the determination of critical temperatures is also difficult 
in equilibrium because of the CSD. Now we have an alternative approach, the STCD, by 
which the critical temperature can be also extracted from the power-law scaling behavior in 
the critical regime f^. Therefore, we can perform comprehensive investigations on statistical 
models by the STCD approach which are independent of the study in equilibrium. 

In this paper we present our numerical study for the g-state Potts model on a 2d triangular 
lattice by using the STCD method. Generally, in short-time MC simulations, there are two 
different schemes to calculate the critical exponent: one is based on the power-law behavior 
and the other on the finite size scaling collapse for different lattice sizes within the non- 
equilibrim relaxation processes. Our attention is specially paid to the universal short-time 
evolution from the fully ordered initial state. The calculations show evidence that there 
exists universal scaling already in the short-time regime where the power-law behavior of 
magnetization and auto-correlation is observed, as well as the finite size scaling collapse 
behavior. These results are then applied to estimate the critical exponents and check the 
corresponding results for different schemes. Our numerical results confirm explicitly the 



universality proposal in the short-time regime ||I6 

2 Model and STCD Scaling 

The Hamiltonian for the g-state Potts model with ferromagnetic coupling (J > 0) defined 
on a 2d triangular lattice is given by 

J N 6 

^= -^EE'^-n-;,, f^M = l, (1) 
^ 1=1 l,=l 

where J2i represents the sum over the lattice sites on a triangular lattice N = L x L, and 
li denotes the nearest-neighbors of site i, whose coordinate number is six on the lattice, cxj 
is a spin variable on every lattice site i. It is known that in equilibrium the Potts model is 
exactly solvable, and the critical points are Jc = 0.5493061 ■ ■ ■ and Jc = 0.6309447- ■ ■ for 
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q = 2 and q = 3 respectively 0. In this work we consider both the q = 2 and g = 3 cases 
and calculate their critical exponents {6, z, f3, u) to investigate the universality explicitly. 

As our starting point we consider an 0{n) vector model [n = 1 for the Ising model) 
with the dynamics of model A and let it suddenly quench from a very high temperature 



with small initial magnetization mo to the critical temperature T^. Janssen, Schaub and 
Schmittmann |^ argued that at the critical temperature dynamic scaling behavior already 
emerges within the short-time regime as, 

M^''\t,T,L,mo) = r'^^/^M(^)(r"t,6^/''r,riL,6^0mo), (2) 

where M*^'^^ is fcth moment of the magnetization, r = (T — Tc) /Tc is the reduced temperature, 
(3 and v are the well known static critical exponents and 6 is a scaling factor, while xq, a new 
independent exponent, is the scaling dimension of the initial magnetization uiq. Here a MC 
sweep over all sites on the lattice is defined as the unit of MC time t. 

For a sufficiently large lattice {L oo) and setting r = 0, 6 = t^^^ in the scaling 
form Eq.(0), the power-law behavior of time evolution of the magnetization at the critical 
temperature can be deduced, 

M{t) ~ mQt\ (3) 

where ^ is a new dynamic exponent which characterizes the universality in the short-time 
regime and is related to Xq by 6' = {xq—(3 /v)/ z. The relation shows that, after the microscopic 
time tmic; the magnetization undergoes an initial increase at the critical point and we can 
easily obtain the result of the exponent 9 based on this power-law form. 

On the other hand, one can also focus on dynamic processes starting from an ordered 
initial state to estimate exponents (3/v and z. Although no detailed analytic study has been 
made about this process, MC simulations in a variety of statistical spin models have shown 



that there also exists a similar scaling relation |18|, |19 



M(^) (t , r, L) = h-^P/^M'^^^ {b-H, b'^/'^T, b-^L). (4) 

At the exact critical temperature(r = 0) and setting b = t^^^, Eq. (^) leads to a power-law 
decay for the magnetization, 

M{t) ~ mo^-^/"^ (5) 
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where M{t) is defined by 



M{t) = ^<E^^W>' (? = 2), (6) 

i 

and 

M{t) = ^ < E(^^.W,i - ^) >> (^ = 3). (7) 

i 

Here N = Lx L is also the total number of spins defined on the lattices and < • ■ • > denotes 
the average over independent initial configurations and/or random number sequences. 

However, the independent determination of 1/z/ seems slightly more complicated than z 
and jS/u. We should differentiate In M(t,r) with respect to r at the critical point 

dr In M{t, T)\r=o = t^'^^'^d^, In M(t', r') |,'=o. (8) 

The exponent 1/z/ can be determined from this power-law behavior if the value of z is 
known. In the practice of our numerical calculations, the differential to r is substituted with 
a reasonable small difference Ar. 

Furthermore the dynamic exponent z can be determinated independently. In order to do 
so, we introduce a Binder cumulant 

M(2) it) 

Here the second moment of the magnetization M^^\t,L) is defined as 

M'^'Ht) = ^<(E^^W)'>> (? = 2), (10) 

i 

and 

M^'Ht) = -^,<{T.i^..m-l)r>, (g = 3), (11) 

i 

respectively. For sufficiently large lattices, finite size scaling analysis shows a power-law 
increase of the Binder cumulant, 

U{t,L) r~.t^/'. (12) 

Based on this feature the dynamic exponent z can be determined and then be used to 
calculate the exponent 1/z/ through Eq.(H). 
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As done by many authors 0, [T^ , the critical exponents can be measured both from 
disordered initial states and ordered initial states. We perform our investigations involving 
these two different initial states for MC measurements to obtain more reliable results. Par- 
ticularly we paid much attention to the scaling forms which describe the dynamic processes 
from the ordered initial state instead of the disordered initial states to avoid the finite mo 
effect on the measured results. There is also an advantage to reduce the statistical fluctua- 
tions for determination of critical exponents and critical temperatures from ordered initial 
dynamic processes. 

Then, there is another scheme to estimate these critical exponents, (3/v and z, based 
on the finite size scaling collapse 0. The results estimated should be consistent with those 
obtained by the power-law behavior described above. This scheme gives an independent 
check whether our dynamic MC simulations are reliable or not. Set r = and mg = 0, 
according to the dynamic scaling behavior Eq. (0), the simple scaling relation for the second 
moment magnetization M*^^^ (t, L) is easily obtained 

M^^\t,L) = b-^^/''M^^\b-H,b-'L). (13) 

Using this scaling relation Eq.(]TB|) for a pair of curves of M^'^\t,L) and M^^\t',L') with 
t' = and L' = L/2, we can estimate the critical exponents /5/z/ and z by searching for 
the best fit of collapse of the curve M^'^\t', L') through a global scaling factor 2'^^/'^ . Also if 
we input the value of z got from the power-law the Eq. (|l^), the value of P/u can also be 
estimated by the scaling relation Eq . ([I3D . 

3 Simulation and Results 



It has been demonstrated in previous works p, |2y, ^ that the results obtained by the 
Metropolis and Heat-bath algorithms are consistent with each other and the latter is some- 
what more efficient than the former. Therefore, in the present paper, the MC simulation 
is only performed by the Heat-bath algorithm. Samples for average are taken over 150,000 
independent initial configurations on the triangular lattices with the periodic boundary 
conditions using L = 32, 64 and 128. Statistical errors are simply estimated by performing 
three groups of averages selecting different random seeds for the initial configurations. Our 
simulation is performed at or near the critical temperature. 

We begin the study by considering the dynamic process starting from random initial states 
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with small mo to determine the critical exponent 6. These disordered initial configurations 
with given values of mo are prepared by the sharp preparation method , and the evolutions 
of the magnetization M{t) are measured as a function of MC time t for mo = 0.04, 0.02 and 
0.01 on the = 128^ lattice. The curves of M{t) in a log-log scale are shown in Figs. 1 and 
2 for the q = 2 and q = 3 model, respectively. Obviously, there exist very nice power-law 
increases of M{t) after tmic — 10 and all the curves are almost parallel to each other. Thus 
the 6 can be estimated from the slopes of the curves in the regime of t=[10, 200]. In Table 1, 
the values of ^ as a function of initial magnetizations mo are presented for q = 2 and q = 3. 
Then the final results of 6, after an extrapolation to mo — > 0, are summarized in Table 2 
both for the models on the triangular and square lattice. It is found that the values of 6 
on the triangular lattice are the same as those on the square lattice. This calculation shows 
first evidence of the dynamic universality in the short-time regime. 

Secondly we study the evolution of magnetizations in the initial stage of the dynamic 
relaxation starting from fully ordered initial states. In Figs. 3 and 4, the power-law behavior 
of the Binder cumulant U{t) for q=2 and q=3 on a lattice with = 128^ is displayed on 
a log-log scale, and the nice power-law of U{t) is clearly exhibited. The exponent z can be 
estimated from the slope of the curves: z = 2.145(3) for q=2 and z = 2.148(4) for q=3. In 
Table 2, the values of z measured for the two models are presented. For a comparison, we 
also list in Table 2 the z values of the corresponding results on the square lattice. Then, 
we investigate the short-time behavior of M{t) starting from the fully ordered initial state. 
The time evolution of M{t) for the two models is displayed in Fig. 5 and Fig. 6 respectively. 
Here it is shown that the finite size effect is small, since there is almost no difference for 
the curves of the different lattice sizes. From the slopes of these power-law decay curves the 
values of the index f3/i>z is determined. Then the values oi (3/v are estimated by input of the 
z values. The results are listed in Table 2 for both q=2 and q=3. Our measured results of 
are more accurate than those obtained in the previous works, compared with the exact 
values of (3/v = 1/8 for the q=2 (Ising) and (3/v = 2/15 for the q=3 Potts model. 

Next the MC simulations for dr InM(t) are carried out by taking the difference At =0.02 
on L=128, 64 and 32 lattices. In Figs. 7 and 8 we plot the evolution of c}T-lnM(t) for the 
two models respectively. The nearly complete overlap of these curves on the different lattice 
sizes indicates that the finite size effect can be ignored. We also notice that the power- 
law behavior is not yet seen before t ~ 20, as the effect of microscopic time scale tmic- 
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Table 1: The measured values of 6 versus the initial mo for the Ising model {q = 2) and Potts 
model (g' = 3) on a = 128^ lattice. The last column gives the results of 6 after extrapolation to 
niQ = 0. 





0.04 


0.02 


0.01 


0.00 


q = 2 


0.183(1) 


0.185(1) 


0.189(2) 


0.191(2) 


q = 3 


0.101(1) 


0.089(1) 


0.082(1) 


0.076(2) 



Therefore, we measure the exponent l/uz from the interval [50,500], and the measurements 
yield the results l/z/=l. 027(6) for the q=2 model and l/z/=l. 223(5) for q=3. These results 
are included in Table 2 and compared with those on the square lattice. 

Finally, we study the time evolution of M*^^-* (t, L) with completely random initial states 
(mo=0) on lattices up to L = 128 to estimate the critical exponent jS/u and z by the finite 
size scaling collapse scheme, Eq.(|l3D, as an alternative way. Figs. 9 and 10 show the time 
evolution of the second moment of the mag netization M^'^\t,L) for L = 128,64,32,16 by 
solid lines. The dots fitted to the lines show the data collapse for the rescaled variables 
L' = L/2, t' = and global scaling factor 2^^/*^, with the standard least square fitting 
algorithm. In other words, we rescaled the M*^^^ (t, L) lines on a larger lattice to a smaller 
lattice L' = L/2 by the scaling factor 2^ and the 2^^/'^ for the pair of different lattice sizes. 
The average values for 2(3 /u and 2; of g = 2 and q = 3 are also given in the Table 2. This 
step gives a further independent check that the STCD method is efficient and the results of 
our dynamic MC simulations are reliable. 

4 Summary and Conclusion 

We have systematically investigated the short-time critical dynamics of the q = 2 (Ising) and 
q = 3 Potts models on the 2d triangular lattices by a large-scale dynamic MC simulation. 
Firstly, by observing the power-law increase of the magnetization M{t) from random initial 
states, the new dynamic exponent 6 has been calculated and its values are the same as those 
for the corresponding models on the square lattice. Then we paid our attention particularly 
to the time evolutions of the Binder cumulant U{t), magnetization M{t) and the derivative 
dr In M{t) from the fully ordered initial state to determine the dynamic exponent z and the 
static critical exponents (3 and u. Finally we checked our results independently by using the 
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Table 2: Results for the critical exponents 9, 2(3/v, 1/v and z for q=2 and g=3 on 2d triangular 
lattice. For comparison, the values on the square lattice are also listedQ. 





e 


2B/u 


l/u 

1 


z 


g=2(Triangular) 










power law results 


0.191(2) 


0.250(5) 


1.027(6) 


2.145(3) 


scaling collapse results 




0.252(2) 




2.153(2) 


(Square, Ref.^) 


0.191(1) 


0.240(15) 


1.03(2) 


2.155(3) 


exact 




1/4 


1 




g=3(Triangular) 










power law results 


0.076(2) 


0.256(6) 


1.223(8) 


2.148(4) 


scaling collapse results 




0.266(2) 




2.191(1) 


(Square, Ref.§) 


0.075(3) 


0.269(7) 


1.24(3) 


2.196(8) 


exact 




4/15 


1.2 





finite-size scaling collapse scheme to the time evolutions of M'^'^\t,L) from the completely 
random initial states. The advantage of the STCD method, by our experience, is that the 
CSD is eliminated because the measurements are carried out at the beginning of the evolution 
process, rather than in the equilibrium where all the correlation lengths are divergent. 

In conclusion, the dynamic relaxation in the short-time regime for the q = 2 and q = 3 
Potts models is studied on the 2d triangular lattice, and it is confirmed that, by comparing 
our estimated results 9, j3 and z/ to those on the 2d square lattice, they belong to the same 
universality class for the corresponding models on different lattices respectively. Our critical 
scaling analysis strongly suggests that the simulation for the dynamic relaxations can be 
used to confirm numerically the universality. Further application of STCD is interesting to 
investigate the spin systems with quenched disorder and quantum phase transitions for 
their universal behavior. 
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Figure 1: Plot of M{t) versus t on a log-log scale for the Ising model on a = 128^ 
lattice. The lines present the different initial magnetizations mo=0. 04, 0.02, 0.01 (from top to 
bottom) . 
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Figure 2: The same as Fig.l, but for the ^=3 Potts model. 
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Figure 3: Time evolution of Binder cumulant U {t) on a log-log scale for the Ising model on 
a = 128^ lattice. 
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Figure 4: The same as Fig.3, but for the q—3 Potts model. 
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Figure 5: Power-law decay of M{t) versus i on a log-log scale for the Ising model with 
mo — 1. The hues for different lattices of L=128,64,32 almost overlap. 
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Figure 7: Curves of 8^- InM(t) on a log-log scale for the Ising model with mo = 1 for different 
lattices of L=128,64,32. It is obvious that the finite size effect seems small, and the power-law 
behavior exists after tmic — 20. 
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Figure 9: Curves of the second moment of the magnetization M^'^'>{t, L) for the Ising model 
with mo = 0. The sohd hues present M^'^\t, L) for different lattices of L=16,32,64,128 (from 
top to bottom) The dotted lines present the rescaled M^^^ {t' , V) . (O denote rescaling 32 — > 
16, + denote 64 ^ 32 and □ denote 128 ^ 64) 




20 40 60 80 100 120 140 160 180 200 




14 



